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We present a method for calculating the transmission spectra, dispersion, and time delay charac- 
teristics of optical-waveguide gratings based on Green's functions and Dyson's equation. Starting 
from the wave equation for transverse electric modes we show that the method can solve exactly 
both the problems of coupling of counter-propagating waves (Bragg gratings) and co-propagating 
waves (long-period gratings). In both cases the method applies for gratings with arbitrary dielec- 
tric modulation, including all kinds of chirp and apodisation and possibly also imperfections in the 
dielectric modulation profile of the grating. Numerically, the method scales as O(N) where N is 
the number of points used to discretize the grating along the propagation axis. We consider optical 
fiber gratings although the method applies to all ID optical waveguide gratings including high-index 
contrast gratings and ID photonic crystals. 

PACS numbers: 42.79.Dj, 42.25.-P, 42.25.-Bs 
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I. INTRODUCTION 

Spectral properties of optical fiber waveguide grat- 
ings 0, Q have typically been addressed with coupled- 
mode theory (CMT) 0, H H @. The theory relies on 
entities with a clear physical interpretation and despite 
its approximations (slowly varying amplitude, small in- 
dex modulation, and synchronous approximations) it has 
proved successful in accounting for the behavior of long- 
period gratings (LPGs) and Bragg gratings (BGs) with 
uniform dielectric modulation, i.e. a dielectric modula- 
tion of the grating with constant period and peak value. 
Over the years corrections to CMT for non-uniform grat- 
ings as well as more general solutions have been addressed 
in a number of papers including Refs. M. M. Wi fToL ITU Il2l 

Others have addressed the problem of gratings with 
arbitrary dielectric modulations through iterative solu- 
tion of integral equations in inverse scattering methods, 
such as the Gel'fand-Levitan-Marchenko coupled equa- 
tions pH , Il4| . The iterative solution of the integral equa- 
tions scales as 0(N 3 ), where N is the number of points 
along the propagation axis, and does not always con- 
verge for strong Bragg gratings. The convergence prob- 
lem can be remedied by using the approximate layer- 
peeling algorithm Il5| which is based on a combination 
of transfer matrices | 16| and a differential scattering for- 
mulation |17| . Originating in geophysics the layer-peeling 
algorithm divides a non uniform grating into a number 
of uniform segments. Numerically the method scales as 
0(N 2 ). 

In this paper we offer a Green's function method 
(GFM) based on the solution of the Dyson's integral 
equation that is exact for coupling of co-propagating 
transverse electric waves as well as counter-propagating 



transverse electric waves. The numerical implementation 
scales as O(N). Our calculation of the Green's function 
formally relies on methods from quantum mechanics and 
the many-body theory of condensed matter systems, see 
e.g. Refs. |la HH Interpreting the Green's function as a 
space(-time) propagator immediately provides us a link 
to the transmission and reflection properties of finite- 
length gratings and thus forms a physical interpretation 
of the Green's function. 

The GFM thus offer an exact integral method, which 
also may be formulated in terms of transfer matrices. 
Furthermore, perturbative calculations may be readily 
be carried out with the Dyson's equation. To demon- 
strate the method we calculate spectra, time delay, and 
dispersion characteristics for uniform, as well as apodised 
and chirped fiber gratings. 

The Green's functions and the Dyson's equation have 
previously been used in electromagnetical systems for a 
variety of problems (see Ref. |2fJ and references therein), 
among those finite size 2D photonic crystals pll 12^ . |23| . 

The paper is organised as following. In Sec. [H] we 
describe the problem of optical waveguide gratings. In 
Sec. lIIII we describe the problem in terms Maxwells equa- 
tions. In Sees. HVfVl we introduce the Green's function 
formalism and the Dyson's equation. In Sec. I VII we ap- 
ply the theory to long-period gratings and in Sec. IV 111 
to Bragg gratings. In Sees. IVIIll and llXl discussions and 
conclusions are provided. 



II. OPTICAL WAVEGUIDE GRATINGS 

The waveguide grating is a scattering region intro- 
duced in the waveguide and is effectively equivalent to 
a one dimensional finite size photonic crystal. Waveg- 
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uides support a number of guided modes; one or two 
fundamental modes and an infinite number of higher or- 
der modes. All modes have a forward and a backward 
propagating part. In the so-called long-period grating dif- 
ferent forward propagating modes couple at a resonant 
frequency. Forward scattering involves a very small mo- 
mentum transfer 8(3 — kAn cS where k = uo/c — 2ir/\ is 
the free-space wavenumber and An cS is the difference in 
effective indices of the resonant modes involved. Since 
(5/3A<3 ~ 1 at resonance, the period A<3 of such a grat- 
ing is necessarily long, hence its name. In Bragg grat- 
ings the forward and backward propagating parts of the 
same mode couple involving a large momentum transfer 
8(3 = 2(3 = 2kn cS and thus Bragg gratings necessarily 
have a short grating period Aq. 

The dielectric modulation seen in Fig. ^ can be per- 
fectly periodic with Aq (uniform grating) or can be quasi- 
periodic (non- uniform grating) . It is possible to alter the 
spectral characteristics of a grating by using apodisation 
and chirp in the dielectric modulation Ref. |24|. Such 
techniques have made grat ings widely used as highly se- 
lective bandpass filters 25] , dispersion compensators and 
pulse compressors [2(| l27j and optical sensors |2^| . With 
apodisation the magnitude of the peak dielectric modu- 
lation varies across the length of the grating. 

With chirp the period of the grating varies across the 
grating. Apodisation and chirp, however, may also ap- 
pear as an unwanted error in an intended uniform grating 
due to fabrication imperfections. 

Side-lobe suppression can be done by an apodisation 
known as raised-Gaussian apodisation |29j | where the di- 
electric modulation has a Gaussian envelope. 

In a linear chirp the grating period varies as Ag(z) = 
Aq(1 + where Cj is the chirp parameter and L is 

the length of the grating. Apodisation can be visualised 
as segments of uniform gratings that are concatenated. 
The spectra are different due to different peak dielectric 
modulation, but share the same resonance wavelength. A 
chirped grating can be visualised as the concatenation of 
uniform grating segments with different grating periods. 



III. GENERAL FORMALISM 

In the following we consider temporal harmonic solu- 
tions to Maxwell's equations with electrical fields of the 
form E(r, t) = E w (r)e~ Iwt , with the subscript to indi- 
cating the frequency dependence. The electrical field is a 
solution to a vectorial wave equation, which has the form 
of a generalised eigenvalue problem [2(J, y(J , 

VxVxE w (r)-fc 2 e r (r ± )E tJ (r)-fc 2 £ G (r)E w (r) = 0, (1) 

where the dielectric functions, e, are 3x3 dielectric ten- 
sors. rj_ is the coordinate vector orthogonal to the axis 
of the fiber which we define to be parallel with the z-axis. 
The subscript V refers to what we will call the reference 
system, which is the bare fiber without any grating. e r 



is thus translational invariant along the z-axis. The sub- 
script 'G' refers to the changes in the dielectric function 
induced by the grating. The grating begins at Zi and 
ends at Zf and £G( r ) is thus zero outside this interval. 
Usually £g( v ) is denoted as 5ec(r) because it is assumed 
that it is small, but here we use the former notation to 
emphasise that the presented theory is formally valid for 
any dielectric profile including high-index contrast pro- 
files since Eq. JTJ is completely general. In the following 
£(j(r) is referred to as the dielectric modulation and we 
will for simplicity omit the subscript on E w . 

In order to solve Eq. (JIJ for a general £g we will use a 
Hilbert function space spanned by the monochromatic 
(fixed k) eigensolutions to to the reference system in 
Eq. IT) 

ES?(r)=ESf( ri )e^, (2) 

where (3 m is the propagation constant of the mode m. For 
the inner product we have (E^\E^) = J drE&>t(r) • 
E^>(r), where f is the Hermitian conjugate (trans- 
pose and complex conjugate) of the field vector. Wc 

choose the eigenmode basis {-E^} —1 to be orthonor- 

mal, (E$\e r \EW) = {^f) 2 8 mni where 8 is the Kro- 
necker delta function. This orthonormalisation follows 
from the fact that Eq. is a generalised eigenvalue 
problem. Any linear operator O on a Hilbert space of 
finite dimension can be represented by a matrix acting 
on a state vector as 6 = J2 mn °mn\E^)(E^\, where O 
has matrix elements O mn = J dr E<, l )t (r)O(r)E<; 1) (r). 

For transverse electric modes, dielectric tensors with 
£xz — £zx — £yz = £zy = 0, and it is possible to integrate 
out the transverse degrees of freedom in Eq. to obtain 

d 2 

J2\ E ) + k 2 i r \E) +k 2 s G (z)\E) = 0, (3) 

where \E^ is a linear combination of eigenmodes with 
scalar coefficients, i.e. (state) vector in the eigenmode ba- 
sis {E$ } _ . We have thus transformed the fully vecto- 
rial eigenvalue problem of Eq. into a scalar eigenvalue 
problem given by Eq. J3J| in the Hilbert function space 
spanned by the fully vectorial eigenmodes. 

IV. GREEN'S FUNCTIONS 

In the following we address Eq. © with the aid of 
Green's functions which will in turn provide us with in- 
sight in the transmission and reflection properties of the 
grating. 

A. Solution with Dyson's equation 

Let G (0) be the zero-order Green's function solution to 
the eigenvalue equation, Eq. © with eq = 0. This has 
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the analytical solution 



and the transmission and reflection probabilities become 



'Cmn( z i z ') — $n 



J0 m \z—z'\ 



(4) 



The zero-order Green's function only depends on z' — z 
in agreement with the assumption of translational invari- 
ance of the reference system. The zero-order Green's 
function can conveniently be split into a spatial forward 
propagating and a backward propagation part 



(z - z')iGf >mm (z, z') + 6{z> - z)iG^ mm (z, z>) 



9{z-z'y- 



if3 m (z~z') 



2f3 n 



+ 8(z'-z)- 



-ifl m {z-z') 



(5) 



The subscripts '+' and '-' refer to the forward and 
backward propagating part of the Green's function, re- 
spectively, and 6 is Heaviside's stepfunction. We also 
note that G { + and G_ have the convenient mathemat- 
ical property that they each split in their variables: 

e ±ip m (z-z<) = e ±i/3 m z e ^i0 m z' _ In terms of G (0) the 

Green's function solution to the full grating problem 
Eq. J2Jl is given by the Dyson's equation [20l | 



T 



Zi)\\ (9a) 



Rrrin — |^rnn| — 2 (3 m [3 n \iGrnn ? %i) \ 



(9b) 



The group delays of the transmission and reflection are 
given by the phase of the amplitude coefficients, t mn and 
r mn . The group delay for the scattering n — > m is given 
by 



arg (tmn) , 



(10) 



where arg(s) is the argument of the complex number of 
x. For the group delay of the reflection scattering similar 
arguments hold and the group delay may be obtained by 
substituting t mn with r mn in Eq. i|l(J|) . The dispersion, 
D, of the n — > m scattering is defined as the derivative 
of the group delay with respect to wavelength 



dX 



(11) 



iG(z,z') = iG {0) {z,z') 

fZf 

— I 



i [ dz"iG {0) (z,z")k 2 e G {z")iG{z",z'). (6) 

J Zi 



The integration interval is changed to [zj, Zf] since eq is 
zero outside this interval. 



B. Physical properties of the Green's functions 

The transmission and reflections properties may be de- 
scribed in a scattering basis formulation where ingoing 
and outgoing modes are related to each other by a scat- 
tering matrix 



S = 



r t' 
t r' 



(7) 



Here, t mn is the transmission amplitude for forward scat- 
tering (left-to-right direction) from mode n to mode m 
and r mn is the reflection amplitude for back-scattering 
(left-to-left direction) from mode n to mode m. Like- 
wise, t' and r' are amplitudes for the opposite directions. 
We restrict ourselves to lossless media for which S is 
unitary and symmetric so that t' = —t', r' = r\ and 
1 = r'r + t't. The Green's function is naturally inter- 
preted as a propagator for field amplitudes and the full 
Green's function may be related to the bare Green's func- 
tion via 

e i{P m Zf-P n Zi) 

iG mn (zf,Zi) — t mn — - — ^^=== — (8a) 



2V PmPn 



(8b) 



V. THE INTERMEDIATE GREEN'S FUNCTION 

The dielectric modulation of the grating not only 
causes the different modes to coupled, but also the in- 
dividual modes to change. This phenomenon is also ref- 
ereed to as self-coupling and has the physical interpreta- 
tion that the grating causes the individual propagation 



constant of a mode (/3„ 



,cir 



k) to change slightly. Self- 



coupling is typically of second importance to cross-mode 
coupling in fiber gratings and it can in many instances 
be ignored. However, we can at no cost incorporate it in 
the free forward and backward Green's function propa- 
gators and we refer to these new propagators G± as the 
intermediate Green's functions. By this trick we eventu- 
ally end up with a Dyson's equation for the full Green's 
function with a particular simple structure. To see this 
we split Sg into two Hermitian operators: a self-coupling 
operator (s), which is diagonal in its indices, and a cross 
(x) coupling operator which is purely off-diagonal: 



£ G (z) = £ s (z) +£ X (z) 



.(*) 



£ mn( Z ) — e G,mn{ z ) £ mn( Z ) 



(12a) 
(12b) 
(12c) 



The diagonal, self-coupling component, e s , is often refer- 
eed to as the 'DC' component and the cross coupling com- 
ponent, e x , as the 'AC' component in analogy with elec- 
tric engineering. We first consider the forward Green's 
function. To zero order in the self-coupling the interme- 
diate Green's function equals the free Green's function, 
G ? \z',z) — G l ^_\z l ,z). Iterating Dyson's equation and 
keeping terms to linear order in s s gives the correspond- 
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ing contribution to the intermediate Green's function 
-i f dz"iG ( P(z',z")k?e B (z")iGf (z",z) 

J Z 



iG f(z',z)^- dz"e s (z"). 



2/3 



(13) 



With the operator 



AP(z',z)-(z'-z) = ^J* dz' 



i%z"), (W) 



we may solve Dyson's equation, Eq. ©, straightfor- 
wardly to any order. One might wonder why we chose to 
define a function as in Eq. 1)14(1. This definition has the 
advantages that the left hand side is distributive in z and 
z', A/3(z', z) ■ (z' -z) = A/3(z', z") ■ (z' - z") + Aj3(z" , z) ■ 
(z" — z), and if the equation is divided by z' — z, then it 
can be seen that Af3 mm , is proportional to an average of 
Eg over z' — z. We will use these properties in the follow- 
ing. For the backward propagating intermediate Green's 
function the analysis above may be repeated. Using the 
distributive property of A/3 and introducing 



0(z)z = f3z + A/3(z, Zi) ■ (z - Zi ) 



(15) 



we may write the forward and backward propagating in- 
termediate Green's functions as 

e 0{z)z p -0{z')z' 

6{z-z')iG + {z,z') = 6{z-z') — , (16a) 

6(z' -z)iG-(z,z') = 6{z'-z) — , (16b) 

zp 

which are diagonal matrices in their mode indices. The 
Dyson equation (JSJ may be rewritten in terms of the in- 
termediate Green's functions and the cross coupling op- 
erator 



iG(z',z) = iG(z',z) 



dz" iG(z', z")k 2 e*(z")iG(z", z). (17) 



A. Numerical evaluation of Dyson's equation 

The Dyson's equation (|17fl describes a system of M x M 
coupled equations. Due to the fact that the matrix iG is 
diagonal in its mode indices the different matrix columns 
of iG are not coupled to each other, and the Dyson's 
equation constitute M problems of size M instead of a 
single problem of size M x M. If only a single elec- 
tric mode is incident on the grating (as is often, but 
not always, the case) only a single column of the matrix 
Dyson's equation needs to be solved. For large M and 
N this gives a significant reduction in time and memory 
consumption. 



The matrix Dyson's equation (|17fl mav be solved 
by simple iteration by inserting G on the position of 
G on the right hand side of the equation to obtain 
ELo *G ,a) (z,z') on the left hand side. This in turn 
may be inserted to find X)a=o iG {c " } {z, z'). Such an it- 
eration may be continued until any order. In this work 
we evaluated Dyson's equation with the discrete dipole 
approximation 31]. In the approximation the scatter- 
ing interval [zi,Zf] is discretised into N cells. Each cell 
is assumed small enough that the variations of the elec- 
tric field E(z) and the dielectric modulation Sa(z) are 
negligible within the cell. 

The numerical task of evaluating a single iteration 
with Eq. l|17fl generally scales as iV 2 where N is the 
number of cells on the discretised z grid. However due 
to the distributive property of the forward and back- 
ward propagating intermediate Green's functions shown 
in Eqs. (|16all6b|l . the calculation scales as N. 



VI. COPROPAGATING WAVES - 
LONG-PERIOD GRATINGS 

In a LPG the power of an incident core mode is cou- 
pled into a forward propagating cladding mode which 
has a different, smaller propagation constant, /?, but the 
same frequency uj. The cladding modes are slightly lossy 
causing the mode power of the cladding mode to slowly 
dissipate into the environment of the waveguide. Thus in 
practice the total energy transmitted through the waveg- 
uide is reduced by the amount of energy transmitted into 
the cladding mode. Physically, the cladding mode losses 
are small over the length of a grating and may be ne- 
glected in the calculations. 



A. Resonance condition 

To analyse the resonance, we may neglect higher order 
contributions than first order in e x (but retain all orders 
in e s ) to the total Green's function 

iG(z f , = iG m (z f , Zi ) + iG w (z f , z^ 

e -i(3(z f )z f ( ^ e 0(zi)z; 

1 - in(zf, Zi)(zf - z^ 



2/3 



where the cross-coupling function, k, is defined as 



2/3 ' 
(18) 



Kmn(,Z , z) ■ (z z) 



„2 r z 



2v f3 m {3 n 



dz 



« e -i(^(*'0-M*'0)*''^ n(2 "). 



(19) 



We can thus estimate the resonant wavelength by finding 
the maximum of |K mn | with respect to the wavelength. 
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For uniform gratings with constant period and strength 
we may rearrange the integrand of the expression as |37j| 

p-i ( An {Z)-Pm (*) -P n (z)+P n (zj) X X ( S 
° c mn\ 4 'J m 

(20) 

The middle factor on the right hand side of Eq. i|20|) 
does not have a different period than e^ n (z). The 
self-coupling thus changes the resonant frequency from 
8/3 m - p n = to an average (3 m (z) - (3 n (z) = |^ in 

the limit ^ 2 ' — > oo. This result is the same as that of 

G ■— 1 1— I 
the synchronous approximation used in CMT 5, 6] which 

arises from different arguments. 



B. Numerical solution of Dyson's equation 



We consider a system of a pair of coupling modes: a 
core mode and a cladding mode. Only the forward prop- 
agating Green's functions need to be considered as the 
coupling to the backward propagating modes is far from 
resonance and is negliable. The use of only forward prop- 
agating Green's functions alters the integration limits in 
the Dyson's equation, Eq. Q17[l. due to the step function 
in the forward propagating intermediate Green's func- 
tion, Eq. H16a(l . We consider a system where only a single 
core mode is incident on the grating and thus need only 
to solve the first column of the matrix Dyson's equation, 
Eq. O 



o 



ik 2 dz' 



iG co {z,z') 
iG ci (z,z') 





"cLcoV^ > 



co,cl( Z ') 





iGc\ x co{% ; z i) 



(21) 



where we have set the notation as G co (z, Zi) = Gn(z, Zi) 
and G c \(z,Zi) = Gii{z, Zi). The contributions to the to- 
tal Green's function from zero to second order in e x can 
be seen in Fig. [5] Numerically, we can improve the dis- 
crete dipole approximation employed in this paper by 
rescaling the top and bottom row by \/G co (z,zi) and 
1/G c i(z, Zi), respectively, since the rescaled total Green's 
function varies more slowly within the cells than the orig- 
inal Green's function due to Aq 3> A for LPGs. To 
solve the system we use the simple iteration method dis- 
cussed in Sec. IV Al The iteration may be continued until 
the desired accuracy is reached. In this work we iter- 
ate until the root mean square of the residue between 
two consecutive iterations is less than 10~ 12 . The di- 
electric modulation is chosen as a uniform modulation 
£g(i") = Ae(r^)i(l — cos(^-)z). We model a typical 
fiber LPG with resonance tentatively at A ~ 1550nm 
with the peak dielectric change Ae = 2.5 x 10~ 3 , a 
grating period of Aq = 500//m, and a difference in the 
effective indices of - nf = 15 %° nm . The differ- 
ence in the effective indices is assumed constant with 
respect to the wavelength. The elements of the opera- 
tor i, are given by f^™W = g gQ fLsiM _ Q 2 Q, and 

e x (r) £ x (r) 

co ' c ,' , = cl|C ? , = 0.10. Finally, the dieletric constant 
of silica is e s i = 2.10. 

With the given parameters we model the spectra of the 
grating with four different lengths L = z£- — 8mm. = 
13.5mm, — 20mm, ^ = 40mm corresponding to 16, 
27, 40, and 80 grating periods, respectively, and k is given 
by the maximum of Eq. HIH|I with respect to wavelength. 

The LPG couples the core mode power into the 
cladding mode with cosine-like amplitude behavior as 



seen in Fig. In passing we note that the LPG ef- 
fectively corresponds to a rotation of the state vector in 
the Hilbert space spanned by {E co , Eci}. The ? rotation 
seen in the bottom left figure in Fig. 0] can be made into 
a 7r rotation seen in the bottom right figure in Fig. 0] by 
doubling the length of the grating. The resonant wave- 
length found by Eq. I|19|l is I6I5.7 nm. This includes the 
self-coupling of the modes which causes a considerable 
shift of the resonant wavelength of 65.7 nm. The time 
delay and the dispersion of the core transmission can be 
calculated with Eqs. HlOlllfl an d is seen in Figs. 03 
The LPG has a low group delay and dispersion making 
it suitable as a dispersion-free wide-band rejection filter. 

If the dielectric modulation has a Gaussian envelope 
the grating is said to be raised Gaussian apodised |29| . 
The so-called taper parameter which determines the 
strength of the apodisation is g = \ . In Fig. |7| is shown 
the spectrum for a grating that is 40 mm long with 80 
grating periods. The apodisation is seen to reduce the 
side-lobes seen in Fig. 0] We also model a linearly chirped 
grating with the chirp constant chosen as c\ = 0.005. 
Thus the grating period varies 0.5% across the length. 
The number of grating periods is 40 corresponding to a 
total grating length of 20 mm and apart from the chirp 
the resulting spectrum would be identical to the one seen 
bottom left in Fig. 0] The spectrum of the chirped grat- 
ing can be seen in Fig. |SJ The 0.5% variation in the 
grating period is seen to have a clear effect on the spec- 
trum. The side-lobes of the same magnitude as those 
of the uniform grating and the resonant wavelength is 
shifted to 1631. 6 nm from 1615. 7 nm even though the 
average grating period has remained the same. 
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VII. COUNTER-PROPAGATING WAVES - 
BRAGG GRATINGS 

As explained in Sec. [H] BGs couple the counter- 
propagting part of the same mode. The coupling of the 
copropagating modes is so far from resonance that it is 
neglected. We will in the following suppress the mode 
subscript. 

A. Resonance condition 

Following the same line of arguments as in Sec. IVI Al 
we arrive at 

«(/,*). fdz" e-^-'V'eV'). (22) 



This expression resembles Eq. <|19ll but has different in- 
tegration limits due to the step function in Eq. I)16bjl . 
For a uniform grating we may again conclude that self- 
coupling changes the resonance condition from 2/3 = j 2 ^ 

to an average 2(3(z) — in the limit -^r 1 — > oo, which 
also is the result of the synchronous approximation used 
inCMT@@. 



B. Numerical solution of Dyson's equation 

We will solve the Dyson's equation (|17|l by direct nu- 
merical solution. Since we consider a system where elec- 
tric mode is incident from the left, we only have to eval- 
uate the first column of the matrix Dyson's equation 



iG- t+ (z,Zi) 



iG+(z,Zi) 




ik 2 



f*,dz'iG + {z,z') 

II s dz'iG-(z,z') 



£+_«) 
et + (z') 'o 



iG+,+(z', Zi) 
iG- <+ (z', 



(23) 



where we have abbreviated G+ = G+.+ and G- = (_?_,_. 
The equation may be solved in the discrete dipole ap- 
proximation by simple iteration as discussed in Sec. IV Al 
The expressions obtained by simple iteration from zero 
order through second order can seen in Fig. [5] in a dia- 
grammatic representation. It is seen that a reflected wave 
has undergone an uneven number of interactions whereas 
the transmitted wave has undergone an even number of 
interactions. This fact is true for all orders of interac- 
tions. 

The iteration may be continued until the desired ac- 
curacy is reached. In this report we iterate until the 
root mean square of the residue between two consecu- 
tive iterations is less than 10 -12 . When simple iteration 
of Dyson's equation converges this may be met in less 
than 20 iterations. Simple iteration often diverges for 
strong BGs. For uniform gratings it is possible to an- 
alyze the divergence. We will not discuss this analysis 
in detail here, but only assert that simple iteration con- 
verges for uniform gratings when \k\L < ir/2 implying 
T co ,co ^ 0.16.. For stronger and non-uniform gratings 
Dyson's equation \\1'M may be solved as a self-consistency 
problem in iG+ +{z? , Zi). In this paper we have employed 
a generalised minimum residue method for stronger grat- 
ings |32| . In our numerical algorithm the generalised min- 
imum residue method automatically takes over if simple 
iteration of Dyson's equation diverges. 

We choose the parameters of the gratings to closely 
resemble a standard fiber BG with resonance tenta- 
tively at A ~ 1550nm, leading to a grating period of 
Aq = 155 ° nm -| . The peak dielectric modulation is Ae = 



2.5 x 10 3 , and the elements of the dielectric operator e 



» 



0.50. 



are given by — : — 

BGs with four different lengths are shown in Fig. El The 
length is given in the terms of the coupling constant at 
resonance, k, given by Eq. (|22|) . The lengths of the four 
gratings are 1.145mm, 2.290mm, 4.581mm, 9.161mm, for 
kL = 1/2, 1, 2, 4, and with 2143, 4285, 8570, 17140, grat- 
ing periods, respectively. Group delay and dispersion of 
the core transmission can be found using the Eqs. I|10lllfl . 
The dispersion and delay for the transmission coefficient 
for the kL = 1 and kL = 4 are shown in Fig. El and 
Fig. El respectively. 

If we use raised-Gaussian apodisation [2^| with length 
L = 9.161mm and the taper parameter, g = 1/4 we get 
the spectrum in Fig. 1131 Apart from the apodisation the 
grating has the same parameters as the uniform grating 
shown bottom right in Fig. EH The apodisation has re- 
moved the side lobes and riddles seen in spectrum of the 
uniform grating, but has lower reflection (higher trans- 
mission) at resonance. The group delay and dispersion 
are also smoother for the apodised grating as seen in 
Fig.Ul 

Linear chirp has been incorporated in the grating seen 
in Fig. El The chirp parameter is chosen as ci — 5 x 10 -5 
and the length of the grating is 4.581mm. The group 
delay and dispersion of the transmission is seen in Fig. 1161 
It is seen that the chirp affects both the spectrum and 
the group delay and dispersion when comparing with a 
uniform grating. 
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VIII. DISCUSSION 

The numerical results of GFM differ from the widely 
used CMY in two respects. The small oscillation seen in 
Fig. 01 does not occur in CMT, since it results from the 
second order derivative which is neglected in CMT. In 
GFM the oscillation occurs for both LPGs and BGs, but 
vanishes with decreasing coupling strength and averages 
out in general. CMT relies on a Fourier decomposition of 
the dielectric modulation. This is a very good approxi- 
mation when the grating is uniform and has a high num- 
ber of periods. An LPG usually have a modest number 
of periods compared with an BG, and the GFM spectra 
differ from CMT spectra in the side-lobes, while they co- 
incide around the resonance. The GFM could therefore 
be applied for better apodisation of LPGs. BGs have a 
large number of periods, and thus CMT and GFM agree 
well for the uniform and the weakly apodised and chirped 
BGs presented in this report. However, BGs are being 
designed with increasing complex dielectric modulations 
to obtain extraordinary spectral and dispersive charac- 
teristics. The GFM may here be useful in the respect, 
that it does not divide the dielectric modulation into uni- 
form segments, and therefore offers more freedom when 
designing the apodised dielectric modulation. The GFM 
could also be used to study the mechanisms of apodisa- 
tion, which are still not completely understood |35| . 

The examples considered have all been fiber gratings. 
The GFM method is exact for transverse electric modes, 
but some optical fibers do not support these transverse 
modes. However the non-transverse part of a mode usu- 
ally only holds a very small fraction of the total electrical 
field and contributes only negligible to the grating char- 
acteristics. 

Although the presented theory is exact from the trans- 
verse electrical wave equation Eq. the theory is in 
fact only exact within the basis set of the electric fields 
of the reference system (in finite element methods this is 
known as Galerkin error being zero) . The electrical fields 
of the reference system may or may not describe the grat- 
ing well. This is actually a feature in all grating theories 
and should not cause any worries, but some consideration 



should be taken when applying the GFM to high-index 
contrast systems. A discussion of the variational error of 
the basis set can be found in |3(t| . 

The Green's functions and the Dyson's equation have 
here been used to model linear dielectric modulations. 
The theory can be adapted to a non-linear dielectric mod- 
ulations, both local and non-local, for simulating har- 
monic generation, soliton dynamics and other non-linear 
effects. If a non-local dielectric modulation is used then 
powerful perturbative methods can be employed such as 
the Feynman diagrams used in many-body physics. 

In Figs. El [HI h is seen that the group delay for both 
Bragg and long-period gratings is negative near reso- 
nance. This may seem conflicting with casuality, but 
physically it is perfectly possible that the group veloc- 
ity (pulse envelope velocity) is greater than the speed of 
light or even zero [3j| causing negative time delays. How- 
ever, numerical simulations indicate that negative group 
velocities are impossible in passive one-dimensional peri- 
odic media psij . 



IX. CONCLUSION 

A method for calculation of optical-waveguide grat- 
ings' optical characteristics based on Green's functions 
and the Dyson's equation has been presented. We have 
presented accurate spectra for long-period and both uni- 
form and apodised BGs as well as the group delay and 
dispersion characteristics. 

The GFM is exact for transverse electric modes for 
gratings with arbitrary dielectric modulations in both 
the long-period and the BG application. The method 
incorporates self-coupling in an exact analytic manner 
and also gives an exact resonance condition for gratings 
with arbitrary dielectric modulations. The GFM relies 
on iterative solution of integral equations which scales as 
O(N) which may be solved to within any given accuracy. 

The GFM may also be useful in simulating non-linear 
problems or imperfect gratings with random mode scat- 
tering. 
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FIG. 1: An optical fiber waveguide with a grating consisting of a dielectric modulation in the core (grey area) with period 
Aq. The incident electric mode in with propagation constant (3 is partially reflected, (r), and partially transmitted (t) by the 
grating. No radial symmetry is assumed in the presented theory. 






FIG. 2: Diagrammatic representation the zero to second order contributions to the total Green's function Eq. 12 U . An incident 
core mode interacts with a copropagating cladding mode through the dielectric modulation. 
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FIG. 3: Uniform LPG at the resonant wavelength. Top: The transmission of the core mode, T co , as function of the grating 
length. Bottom: the uniform dielectric modulation of the grating. 
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FIG. 4: Spectra of uniform LPGs with a uniform dielectric modulation (bottom in Fig. E| . The spectra is plot calculated for 
four different for grating lengths, Zf — Zi = y^-, jjj, jjjj, ~ corresponding to 16, 27, 40, and 80 grating periods, respectively. At 
kL = 7f (bottom left) all core mode power is transferred into the cladding mode and back again into the core mode at kL — n 
(bottom right). 
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FIG. 6: Time delay and dispersion of the uniform long-period grating with length Zf — Zi — —. The corresponding spectrum 
is seen in the bottom right of Fig. 2] 



12 




10 20 30 



z[nm] y 



1350 1400 1450 1500 1550 1600 1650 1700 1750 1800 1850 1900 

X [nml 



FIG. 7: Apodised LPG. The side-lobes seen in the bottom left of Fig. [I] have been reduced. Insert: a schematic view of the 
raised Gaussian apodised dielectric modulation. The grating period has been exaggerated with respect to the grating length 
for clarity. 
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FIG. 8: Linearly chirped LPG. The side-lobes seen in the bottom left of Fig.[4]have been enhanced rather than reduced. Insert: 
a schematic view of the linearly chirped grating period of the dielectric modulation. The grating periods and the chirp have 
been exaggerated with respect to the grating length for clarity. 
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FIG. 9: Diagrammatic representation of the zero to second order contributions to the total Green's function Eq. 1231 for a BG. 
An incident core mode interacts with its oppositely propagating counterpart through the dielectric modulation. 
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FIG. 10: Uniform BG spectra for four different lengths, nL = 1/2, 1,2,4. The Bragg rejection band is emerging for increasing 
kL. 
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FIG. 11: Group delay and dispersion for the transmission coefficient of a uniform BG with kL = 1. The corresponding spectrum 
is shown top right in Fig. I1UI 
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FIG. 12: Group delay and dispersion for the transmission coefficient of a uniform BG with kL = 4. The corresponding spectrum 
is shown top right in Fig. 1101 
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FIG. 13: Spectrum for a raised- Gaussian apodised BG. Insert: the raised- Gaussian dielectric modulation. The grating period 
has been exaggerated with respect to the length of the grating for clarity. 
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FIG. 14: Group delay and dispersion for a raised-Gaussian apodised BG. The spectrum of the grating is shown in Fig. 1181 
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FIG. 16: Group delay and dispersion for a linearly chirped BG. The spectrum of the grating is shown in Fig. 1151 



